235 


Mastic Solids at Best or in Motion in a Liquid . 

opinion, the view I put forward in 1877 that “ new stars ” are produced 
by the clash of meteor-swarms; and they have suggested some further 
tests of its validity. 

We may hope since observations were made at Harvard and Potsdam 
very near the epoch of maximum brilliancy, that a subsequent complete 
discussion of the results obtained will very largely increase our know¬ 
ledge. The interesting question arises whether we may not regard the 
changes in spectrum as indicating that the very violent intrusion of the 
denser swarm has been followed by its dissipation, and that its passage 
has produced movements in the sparser swarm which may eventuate in 
a subsequent condensation. 

My best thanks are due to those I have named for assistance in this 
inquiry. 


“ Elastic Solids at Eest or in Motion in a Liquid.” By C. Chree, 
Se.D., LL.D., E.R.S. Received November 19,—Read Decem¬ 
ber 13, 1900. 

§ 1. The problems dealt with in the present paper are probably of 
little practical importance ; but they seem of considerable interest 
from the standpoint of dynamical theory. The hard and fast line 
which it is customary to draw between Rigid Dynamics and Elastic 
Solids has been discarded, and a more direct insight is thus obtained 
into the modes of transmission of force in solids. 

Let us consider a solid of any homogeneous elastic material, possessed 
only of such symmetry of shape as will ensure that if it falls under 
gravity in a liquid, each element will move vertically. Take the 
origin of rectangular Cartesian co-ordinates at the centre of gravity, 
the axes of x and y being horizontal, and the axis of z being drawn 
vertically downwards. At time t let £ be the depth of the C.G. below 
a horizontal plane in the liquid, the pressure on which is uniform 
and equal to II. The existence of gaseous pressure on the liquid 
surface would only contribute to n without modifying the general 
conditions of the problem. 

Consider first the elementary hydrostatical theory, according to which 
the liquid pressure at any point x, y,-z on the surface of the solid acts 
along the normal, and is equal to 

n +gp'(C+z), 

where p is the density of the liquid, supposed uniform. 

If the solid fall or rise very slowly, and the viscosity of the liquid 
is very small, the results based on the hydrostatical theory ought to 
give a close approximation to the truth. 
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If a, J3, y represent the elastic displacements, xx, xy, &c., the stresses 
in the notation of Todhunter and Pearson’s ‘ History of Elasticity/ 
the body stress equations are of the type 


dxx 

dxy 

dxz 

d 2 oi 

= o, 

dx 

dy 

+ dz 

~ p dP 

dxy 

d <iy 

dyz 

d 2 f3 


dx + 

dy 

+■ clz 

~ p dt 2 

= o, 

dxz 

dyz 

dzz 

+ P^9- 

d 2 (C+y) 

dx + 

dy 

+ ~dz 

dP 


(i); 


where p represents the density of the solid, g the acceleration of 
gravity. 

The equations treat x, ?/, z as constants for each element of the 
solid, and so assume that the motion, if motion takes place, is 
purely translational. 

If X, /x, v be the direction cosines of the outwardly directed normal 
at a point x, y, z, the surface equations are 


(Xxx + gxy -f vxz) / X = (Xxy -f gyy + vyz)/g = (Xxz + gyz + vzz)/v 

= H -9P'(C +Z ) :. 


( 2 ). 


The equations ( 2 ) are satisfied by the assumption 
S' = yy = zz - -IL-gp'(£+z), 
xy = xz = yz = 0 


(3), 


Also the values (3) satisfy the body stress equations (1), provided 


d?cc\_ <P/3 
dP ” dP 


y) 

'l w 




= ~9P 


(4). 


We can satisfy (4) by assuming 



C = const. + ^<7 ^t ft 
P 


(5). 


For brevity, the constant in (5) will be supposed to be zero. 

The result (5) is of course that given by ordinary elementary 
methods for the accelerated motion of a solid rising or falling in a 
liquid of different density. 
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On looking more closely into the matter an inconsistency manifests 
itself. Supposing for mathematical simplicity that the solid is 
isotropic, of bulk modulus Jc, we find that the displacements answer¬ 
ing to (3) are given by 

ctjx = /% = ~{H + gp'(C+si)}j3Jc, 

y = -[IIz+$'/>'{*(£+«)- J(® 2 + y 2 + « 2 )}]/3i 

The inconsistency consists in the fact that, by (6), a, (3 , y contain 
terms in £, and so by (5) terms in t 2 , while above it was assumed 
that d 2 a/^ 2 , &c., vanished. It thus appears that the solution embodied 
in (3) and (6) is valid and complete only when £ does not vary as t 2 , 
i.e., only when the solid is at rest or moving with uniform velocity in 
the liquid. 

Though thus restricted, the solution is notable from its simplicity 
and generality, as applicable to any homogeneous solid (free from 
cavities) at rest in a liquid of equal density. 

The values (3) for the stresses apply irrespective of the species of 
elasticity. The displacements are given by (6) only when the material 
is isotropic, but corresponding expressions are immediately obtainable 
for materials of greater complexity. If for instance we have material 
symmetrical with respect to the co-ordinate planes, we have 

-x{Il + gp(£+z)}(l - ^ 12 ~i?i3)/Ei, 

P = - y {n+#>'(£+s)} ( l-ysi- V2s)/Es , 

... (7). 

y .= - * {n + gp (C + i*)}(l - mi - W/Es 

+ - Vu - Viz) + 1 ^- (! - V21 - Vis) ) 

Here Ei, E 2 , E s are the three principal Young’s moduli, while 
r /i 2 j ^ 13 ) &c., are the corresponding Poisson’s ratios. 

§ 2. Presently we shall consider the equilibrium problem in greater 
detail. Meanwhile, in the case of uniformly accelerated motion, we 
shall obtain a self-consistent solution for a sphere, or any form of solid 
ellipsoid, under the conditions assumed in § 1. 

The procedure to be adopted is the same for all species of elastic 
material. If for definiteness we suppose the material symmetrical with 
respect to the three co-ordinate planes, we first assume that the 
stresses (3) and displacements (7) form part—but only part—of the 
complete solution, £ being given by (5). Then substituting from (7) 
in the body stress equations (1), we find that the stresses of the 
supplementary solution, as we may call it, must satisfy 
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where 


dxx clxy dxz 
dx dy dz 

dxy dyy dyz 
dx dy + dz 

dxz dyz dzz 
dx + dy + dz 


- Fpx , 

-Q py> " 

- Hpz 


P p = if-p (p-p ')(l - m2 - ms)l^i, ) 


Qp — g 2 P '(p -p')(l- V21 - V2a )/E. 2 , • 
P P = fp (p~p)(l-~Vs l-^32 )/Es . 


(8); 


(9). 


The surface equations to be satisfied by the supplementary solution 
are 

Xxx + yo:y 4- vxz = A#?/ 4- pyjy 4- vyz — te? 4- pyz 4- vzz = 0 ... (10). 

The problem thus resolves itself into that of an ellipsoid acted on 
solely by bodily forces derivable from the potential 

J(Pa 2 + Q/4-R? 2 ). 

This problem was solved by me in 1894 for isotropic* materials, and 
in 1899 I extended the solution to seolotropicf ellipsoids. We can 
thus derive a satisfactory supplementary solution from the sources 
specified. Finally adding the stresses of the supplementary solution 
to the stresses (3), and the displacements to the displacements (7), we 
have a consistent and complete solution of the problem presented by 
a heavy ellipsoid in a homogeneous liquid, when the action of the 
liquid is supposed that given by elementary hydrostatics. 

§ 3 . The supplementary solution, though simple in type, contains 
terms which are of great length when the ellipsoid has three unequal 
axes, and is of a complex kind of seolotropy. It will thus perhaps 
suffice to select for illustration the simple case of an isotropic sphere of 
radius a. 

Denoting Young’s modulus by E, Poisson’s ratio by ?/, and writing 
r 2 for x 2 4 - y 2 4 - z 2 , we have in full 


# ‘Roy. Soc. Proc.,’ yol. 58, p. 39; ‘ Quarterly Journal of Pure and Applied 
Mathematics,’ yoI. 27, p. 338. 
f ‘ Camb. Phil. Soc. Trans.,’ vol. 17, p. 201. 
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£»= -H -gp’(g + yZ—Lf) 

P 

+{(s - -'?) (« 2 -0-0 +s-j) (j/ 2 +o }, 



-Il-r Jp ’(z + lgtULp) 

P 

- o+w+ot 



-n-«>'(*+fc£z£'f*) 

p 

+ 3 - W- O - (1 + 3,) (** + y% 

#?//&$ ~ xz/oxz = yzjyz = - # 2 />'(p - p )(1 - 2??) 2 ~ {5E(1 - *))} 


*/* = m = - I^{n+^v+k^ 2 )} 

y = - Ig- {n^+ (/p'«(s+ P -EJL _ x^',.2j 


( 12 ). 


The terms in g 2 constitute what has been called above the supple¬ 
mentary solution. In the case alike of the stresses and of the dis¬ 
placements they are exactly the same as if the sphere were under a 
self-gravitative force which followed the ordinary gravitational law, 
and which had for its accelerative value at the surface of the sphere 


ffV (P - P) 1 - 2 »> 


This imaginary gravitative action represents attraction or repulsion 
between elements of the solid according as p~ p is negative or posi¬ 
tive. It is thus an attraction when the sphere rises in a heavier liquid, 
a repulsion when it sinks in a lighter. The smaller I - 2?;, or in 
general the less compressible the solid, the smaller is the effect of this 
imaginary gravitative force relative to that of the hydrostatic pressure 
II 4 ~ gp{z 4- {); on the other hand its relative importance increases 
rapidly with the size of the sphere. 

Representing by dashed letters the parts of the displacements 
depending on p ~ p, we have 
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a 'M = P'/y = y’/z 

~ g2p \i l) ~ (1 - l). - V)0?~ (1 + r,y~] - < 2 ] (13). 

At the very beginning of the motion, the expression inside the square 
bracket is positive for all values of r ; but as t increases it changes 
sign, first at the surface, last close to the centre of the sphere. If £ at 
C o represent the distances fallen when the expression vanishes at the 
surface and at the centre respectively, we have 


Ca/a = (l- v )g(p-[>')aim,-\ 

^ .( 14 ). 

Co/« = (3-v)g(p-p)a/30k J 

Unless a is enormously large, ( a and (o must be extremely small for 
any ordinary elastic material. 

In reality, in order to be instantaneously at rest, the sphere would 
require to be supported or acted on by some suddenly suppressed force, 
or to be in the act of reversing some previously impressed motion. 
The elastic strains and stresses might initially retain the impress of 
the pre-existing state of matters, and there are thus special sources of 
uncertainty affecting the applicability of (14) to actual conditions, 
which should not be lost sight of. 

§ 4. The problem just considered has been advanced as showing how. 
under a consistent dynamical system, producing uniform acceleration 
in a straight line, there appear elastic strains and stresses which simu¬ 
late the action of self-gravitation in the material in motion. The 
conditions postulated do not answer exactly to what happens when a 
real solid moves through real liquid at the earth’s surface. Under 
such circumstances the action between solid and liquid is not fully 
represented by the hydrostatic pressure. If the fluid be “perfect,” 
ordinary hydrodynamical theory* gives for the pressure p on the 
surface of the sphere, supposing u the velocity, 

p = II -f gp(C + s) + p'daiiPx + f #P 2 - .(15), 

where Pi, P 2 are zonal harmonics, whose axis is the vertical diameter. 
We shall now consider this case, on the hypothesis that the velocity is 
so small that teiuns in u 2 are negligible. Instead of (3) and (6) we 
find for the stresses and displacements, the material being supposed 
isotropic, 

xx - yy « zz « -ll-gp'(t+z)-iup'z 9 

✓-V V V 

xy = xz = y z ~ 0 



# Cf. Lamb's ‘ Hydrodynamics/ Art. 91. 
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cc/x = Sly = - {II + gi>'(t+z) + 


y = _ 

* IV 


=p[n *+9M+y<9 +i“) (* 2 -y 2 )] j 


{■••• (17). 


Instead of (4) we have 

Also u EE dH/dt*, 

thus, if d 2 y/dt 2 be omitted, we have 


(p + iP) ; 


g{p-p), 


or £ = constant + 4a —— ~t" .(18). 

' p + ip 

This is, of course, only the well-known result, that the dynamical 
action of the liquid may be regarded as adding to the mass of the 
sphere that of a hemisphere of the liquid.* We may suppose the com 
atant in (18) to be zero, suitably interpreting II. 

As in the first case considered, the existence of t 2 in £ and, conse¬ 
quently, in a, yd, y, makes a supplementary solution necessary. The 
stresses of the supplementary solution must satisfy the surface equa¬ 
tions (10) as well as the following body stress equations : 

/ + fay + ( d 3. a. j. <&] /,, 


dx + dy dz JI * \dx + dy + dz j / 


dxz dyz dzz\ / ^ 
dx dy + dz j / 


1 - 2 t? 2 p V(/> -/) 


.(19). 


It will be observed that the retention of the term in u in the pres¬ 
sure has only modified (reduced) the acceleration without altering the 
type of the supplementary solution. It will thus suffice to record the 
complete expressions for the displacements, viz., 




_M_m _,.c . , g 2 p'(p -p ') ?f il 

2(2p + p') ( 2 P + p' 


(1 - 2rjf gW(p ~ />'), 

5(1-^)E 2 2p + p 


[(3-r?)a 3 -(l+ij)r 2 ] 


* fy. Lamb’s ‘ Hydrodynamics/ Art. 91; or Basset’s ‘Treatise on Hydro¬ 
dynamics/ Art. 182. 
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In obtaining this solution we have neglected terms in u 2 } i.e terms 
in (d£jdi) 2 or g 2 t\p - p) 2 j{p + ip) 2 i in the expression (15), while there 
appear in the solution terms containing g 2 t 2 (p — p)f(2p + p). Thus 
our work is consistent only when (p - p)/p is small, and even when 
this is the case the fact that u 2 increases as t 2 involves a restriction 
which should not be overlooked. It would not, I think, be a very 
difficult matter to obtain a complete solution answering to the full 
value (15) of p. Treating u 2 at first as a constant, we could at once 
write down, from my general solution* for the isotropic elastic sphere, 
the displacements answering to the surface pressure |p% 2 ( 3 P 2 — 1): 
but the explicit determination of the corresponding supplementary 
solution would be much more laborious than in the first case treated 
above. 

§ 5. When p and p are equal, and u 2 is thus really constant, the 
complete values of the stresses and displacements answering to the 
surface pressure (15) are as follows :— 


» = - n - gp’(z + ut) + \pu* +1 [(7 + 2 V ) a a 

+ 3rj (x 2 + oy 2 ) - 3 (7 + 6ij)2 2 ], 


7y = - n - 9pip + «0 + Ip u ~ +1[(7 + 2rj) a 2 

+ 3ij (5x 2 + «/ 2 ) - 3 (7 + 6>j) « 2 ], 

« = - n - gp\z + ut) + |p'« 2 - [2(7 + 2 -g)a 2 

-3(7 + 7/) (x 2 + if) + 6r)z*], 

xy ~ ~ 9pit 2 r]xya~ 2 — [2(7 + 

s —v ✓—S 

xzjxz — yzjyz = 9pu 2 ga~~ 2 [4(7 + 5?/)] 


[... ( 21 ); 


ct/x = Ply = ■ - [n - Ip'm 2 + gp\z + ut )] 

+1 ^fT5i~ 2[(7+H ft2 - Cjv {x 2+ f) “ 3 (7 ~ 8r '^> 

7 = -7=^[(U-\p\P)z + gp'{idz+^-x‘>--f)}} 


}• (22). 


- [2 «,) <* + *>-!2,, 2] j 

§ 6. In real liquids viscosity is more or less present, and as the 
hydrodynamical equations have been solved for the case of an ellipsoid 

* i Camb. Phil. Soc. Trans.,’ vol. 14, p. 250. 
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when the retarding action of viscosity neutralises the acceleration due 
to gravity, it is worth considering. The hydrodynamical solution 
really assumes the velocity to be small, and the ellipsoid to be so 
remote from the surface and other boundaries as to be practically in an 
infinite liquid. 

It is not very difficult to deduce from the formulae in Lamb's 
4 Hydrodynamics/ Art. 296,—though I have not seen the result noticed 
—that the viscous surface action reduces to a force fvr per unit surface, 
opposite to the direction of motion, tr being the perpendicular from the 
centre on the tangent plane, and / a constant. The recognition of 
this fact saves us from the labour of considering the general expressions 
for the hydrodynamical pressures, which are of a very complicated 
nature. 

As the motion is steady, the body stress equations are 


/■N /—v -—- -—\ •—X /-N -—x /—>. 

dxx dxu dxz dxu dyy dyz dxz dyz dzz 

tlx + dy ^ dz ” dx + dy + dz ~~ dx dy ^ dz ~ 


0...(23); 


while the surface equations are— a, b , c being the semi-axes of the 
ellipsoid—- 


a~~ 2 xxx + b~~ 2 yxy -f c~ 2 zxz = - a~ 2 x{n + gp'(£+z)}, 
a-3x%y + b-*ijyy + c-*zyz = -b- 2 y{ll + gp({+z)}, 
a~ 2 xxz + b~hjyz + c~ 2 zzz ~ n + gp (f + z) } -/ j 

The surface equations are satisfied by 


(24). 


XX = - n - Cjp'{( + z) + (a-j'?)fz ,' 

yy — -n-g P '(£+z)+{Vi<?)fz, 
zz = - n - gp((+z) -fz, 
xy = 0 , 

xz/z = yz/y = -/ 


(25). 


The values (25) also satisfy the body stress equations (23), pro- 


vided 

~ 3f + 9(p- p) — 0 . 

.(26). 

As 

f j fwd& = 3/. imbc, 



when the integral is taken over the surface of the ellipsoid, (26) is 
simply equivalent to the condition that the motion is not accelerated, 
or that 
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where u is a constant. As to the value of u t it has been proved that 
the total viscous resistance to the motion is* 


l&TT(jfuabc/(xo + c 2 y 0 ), 
where p! is the viscosity, and 

Xo EE abc 1 [(a? -f X) ( b 2 + X) (c 2 + X)]-*dX, 

Jo 

y 0 = J ”[(«* + X) (52 + X) (f + A)*]~yX. 

But this resistance is also equal to g (p - p) tt abc [or to jT/WSJ. 
thus 

u = g{p - p) (xo + c 2 yo)/12/t'. 

Substituting for C and / in (25), we have 


££ = - n - gp(z -f ut) + lg(p - p) a 2 z/c 2 , " 

yy = 

^ =•- -n-g p ’(z + ut)-y(j>-p)z, J-.(27). 

Q = 0, 

= r4 l J = -y(p-p') 

The corresponding displacements, supposing the material isotropic, 


- h^*[n 4-flrp'(* + id)] +■ ® -J + ^ ( 1 “ ?)] ’ 


- ^b[ n + 9pi? +«0]+^ 7 + ’?( 1 


1-977 

1 - " tit 


[II0 + gp utz + \gp{f - a:* - 


| v °L±V 
6E V V 


g(p-~p')r 0 f n „ a 2 ~~7]h 2 \ 0 / 0 5 2 ~r/a 2 \~l 

‘ 6E L*" ( 2 + 3v+ ~ir~J + r (i + otj + — jj 


§ 7. The terms inside the first brackets in (28) contain II or gp , and 
represent displacements which vary only with the depth of the element 
or its distance from the centre of the ellipsoid. The terms containing 

* Cf. Lamb’s c Hydrodynamics/ Art;. 296. 




Elastic Solids at Best or in Motion in a Liquid . 245 

g (p ~ p)i 011 the other hand, depend largely on the shape of the 
ellipsoid. 

Thus, denoting them by a, /3', y', we have approximately, in the 
ease of a very elongated ellipsoid, whose long axis is vertical, 


of/x = Ply = g(p-p') v z/3E, "| 

> .(29); 

y = - (p - p') [« 2 + (2 + 3tj) (® 2 +»/)]/ 6 E J 

and, except in the immediate vicinity of the central section 2 = 0 , we 
ma}^ take in place of (29) 

a'lxr, = p/yrj = -y '/(£») = ?(p-p>/3E.(30). 

In a very flat ellipsoid, approximating to a disc, with the short axis 
vertical, we have approximately 


OL = g(p — p')o; 2 (a 2 ~ ^J 2 )/(3Ec 2 ), 

P = ^(p-p')f(S 2 -V 2 )/(3Ec 2 ), 


S' (31)- 


y' = - g(p — p) [(a 2 — rjb 2 ) o: 2 + (ft 2 - rja 2 ) ?/ 2 + rj ( a 2 -f 6 2 )2 2 ] ~ ( 6 Ec 2 ) 


Except close to the vertical diameter, the terms in z 2 in y' would be 
relatively negligible, while, in general, a! and /3' would be small com¬ 
pared to y'. 

In the case of the sphere it is perhaps more convenient to record the 
complete solution, viz., 


XX = yy = ~Il-gput + lg(p -4p>, 
£ = -Il-^p't^ ~l-g(p + 2p')z, 

xij == 0, 

£/» = yzfy = -lg(p- p ') 


(32); 


*/* = P/y = -i^[n+ ?/ .'(*+irf)]+^(p-p>/B, 

7 = -- : ^p[(n + !;p'M0«+i5'p'(* 2 -a ;3 --y)] 

- [(3 + 2r,) (** + 2/ 2 ) + (1 + 2 V ) * 2 ] 


[March 13, 1901.]—The paper as originally presented to the Society 
dealt briefly with two or three other details. It showed how the solu¬ 
tion in § 6 depended not on the viscous resistance varying as the first 
power of the velocity in the final state, but on its varying over the 
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surface as the perpendicular on the tangent plane. In particular, if, in 
accordance with Mr. Allen’s experiments,* there be possible forms of 
final uniform motion for a sphere in which the resistance varies as 
u% or u 2 (it, being the velocity), it was shown that the solution would 
still be of the form of (32) and (33), provided the distribution of the 
viscous resistance happens to remain unchanged. 

It was pointed out that in an isotropic solid, free of cavities, at rest 
in a liquid, the stresses are everywhere the same as if each element 
were separately subjected to the pressure answering to its depth • but 
that when cavities exist in the solid the state of matters is altered. As 
an example, a complete solution was given for a hollow spherical shell 
fully immersed. 

It was shown that, in a completely solid body, the greatest strain 
and maximum stress-difference theories agreed in indicating no ten¬ 
dency to rupture, but that when cavities existed, it was otherwise; in 
particular, that in the spherical shell there is on either theory a 
tendency to rupture, greatest at the lowest point, which approximately 
in a thin shell varies directly as the depth and inversely as the thick¬ 
ness of the shell. 


“ On the Heat dissipated by a Platinum Surface at High Tempera- 
tures. Part IY.f—High-pressure G-ases.” By J. E. Petavel, 
A.M.I.C.E., A.M.I.E.E., John Harling Eellow of Owens 
College, Manchester. Communicated by Professor Schusteb, 
F.B.S. Beceived February 7,—Bead March 7, 1901. 

(Abstract.) 

The rate of cooling of a hot body in gases at pressures up to one 
atmosphere has received considerable attention, but with regard to 
gases at high pressures practically no data were up to the present 
available. It was thought therefore that an experimental investigation 
of the subject might prove of some interest. 

The experiments were carried out with a horizontal cylindrical 
radiator contained in a strong steel enclosure, the enclosure being 
maintained at about 18° C. by a water circulation. 

It is shown that the rate at which heat is dissipated by the radiator 
may be expressed by the following formula— 

E = ap a + bpP J, 

where E = emissivity in C.G.S. units — total amount of heat dissi- 

# ' Phil. Mag./ September and November, 1900. 
f For Parts I, II and III see ‘ Phil., Trans./ A, vol. 191, p. 501, 1898. 



